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We model the effects of atomic thermal motion on the propagation of a light pulse in an electromag- 
netically induced transparency medium by introducing a set of effectively temperature-dependent 
parameters, including the Rabi frequency of the coupling field, optical density and relaxation rate 
of the ground state coherence, into the governing equations. The validity of this effective theory is 
verified by the close agreement between the theoretical results and the experimental data. 
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I. INTRODUCTION 

The ability of controlling light is crucial to the ful- 
fillment of quantum communications and all their prac- 
tical applications. It has been demonstrated recently 
that a light pulse can be slowed down [IJ [2] and even 
stored in media [31 H] by using the effect of electromag- 
netically induced transparency (EIT), a nonlinear optical 
phenomenon that renders an opaque medium transparent 
by irradiating it with an electromagnetic field [5 j . Subse- 
quent studies indicate that the slow light (SL) effect can 
greatly enhance optical nonlinearity at the low light level 
[5] , and has potential applications in the development 
of low-light-level and single-photon devices, such as the 
all-optical switches [7| and quantum phase gates based 
on the cross-phase modulation scheme [HI [H] • The nov- 
elty of EIT thus provides a new fashion to manipulate 
the behavior of light, and has attracted a great deal of 
research interest. 




FIG. 1: (Color online) The energy levels of a 3-level atom of 
A-type. The transitions |1) «-)• |3) and |2) «-)• |3) are driven by 
laser fields with Rabi frequencies fi p and Q c correspondingly. 



The prototype of EIT systems consists of a weak probe 
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beam of central frequency oj p and a strong coupling beam 
of frequency lu c interacting with a three-level atom, as 
shown in Fig.l. The dynamics of the slowly varying den- 
sity matrix elements in the weak probe limit are described 
by the optical Bloch equation [TU] 
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where f2 p and 57 



denote the Rabi frequencies driving the 
transitions |1) <H- |3) and \2) <H- |3) respectively; 7 is the 
relaxation rate of the ground-state coherence, and T is 
the spontaneous decay rate of the excited state |3). The 



detunings A p and A c 



are defined by A p = UJ31 



and 



A c = w 32 - 



whcre u>ij denotes the transition frequency 



between the energy levels \i) and \ j). The propagation of 
the probe pulse is described by the Maxwell-Schrodingcr 
equation 
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where rj measures the opacity per unit length of the 
medium and is related to the optical density of the 
medium by d opt = 2Lrj/T (L is the medium length). Here 
P31 corresponds to the slowly varying amplitude of the 
optical coherence of the probe transition. 

We note that a decay constant 7 is introduced in eq. 
to account for the relaxation of pi\- Phenomenologi- 
cally, we can write 7 = 7t + 70, where 7t represents 
the effects originating from the thermal random motion 
of atoms and 70 represents other intrinsic effects, such 
as the linewidth and frequency fluctuations of the laser 
fields, the stray magnetic field, and the collisions between 
atoms, etc. Except the demonstration of SL in Bosc- 
Einstcin condensate [HI HI] , most EIT experiments have 
so far been carried out using either laser-cooled atoms or 
hot atoms, with temperatures widely ranging from mil- 
liKelvins to room temperature [131 1 14] . It is recognized 
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that the randomized Doppler shifts stemming from the 
thermal motion of atoms is the predominating mecha- 
nism of decoherence at finite temperatures. Since the 
temperature as an independent variable does not enter 
cqs.Q — |3]), 77- is bound to depend implicitly on the 
temperature and can only be determined by numerical 
fitting with the experimental data. 

In order to gain more insight into the dynamical prop- 
erties of an EIT medium at finite temperatures, it is de- 
sirable to take the atom's motion into account by statis- 
tical approach. Javan et. al. suggested that the macro- 
scopic polarization, which serves as a response to the 
randomized frequency shift, should be obtained by tak- 
ing average over all possible velocities characterized by 
the Maxwcll-Boltzmann distribution [15, 16 . It is ex- 
pected that such an averaging process smears the atomic 
coherence and thus contribute largely to the relaxation 
of pi2- More recently, the dynamics of SL, light stor- 
age and stationary light pulse in finite-temperature EIT 
media are studied in [17j [18] . Wu et. al. [17] have 
developed another method to solve the dynamics of sta- 
tionary light pulses in a cold atomic cloud by solving 
the Maxwell-Liouville equations which contain the higher 
order spin and optical coherence and the Doppler ef- 
fect arising from the atomic thermal motion. By using 
the gauge-invariance of Schrodinger equations under the 
Galilean transformation, Su et. al. [18] have developed a 
numerical scheme to study the dynamics of SL and light 
storage in a Doppler-broadened EIT medium, in which 
the effects of the atom's external degrees of freedom are 
fully considered. 

In this paper, we investigate the dynamics of SL in a 
Doppler-broadened EIT based on the formalism of opti- 
cal Bloch equations. The central idea of this work is to 
treat the atom as a stationary object but otherwise in- 
troduce some temperature-dependent parameters which 
are responsible for the manifestation of thermal effects. 
Specifically, we analyze the output probe pulse under 
the influence of the atomic thermal motion by endow- 
ing the coupling field, the optical density and the phe- 
nomenological decay constant with temperature depen- 
dence, namely, we let tt c fi c (T), 77 77 (T), 7 -> 7 (T) 
for a stationary EIT medium. Derivation of these effec- 
tive parameters will be described in the following. 



II. FORMALISM 

To begin with, we assume that the shape distortion of 
the output probe pulse can be treated as a superposi- 
tion of plane waves of different frequencies propagating 
through the highly dispersive and slightly absorptive EIT 
medium with length L. Each output plane wave attains 
a factor of e l27TL ( n n+ ln i)/ x while propagating inside the 
medium, where ni and hr are the real and imaginary 
parts of the refractive index of the medium, which lead 
to the attenuation and phase shift of the output plane 
waves, respectively. Note that the refractive index is re- 



lated to P31 by 2n {ur + inj) /r/A = p 3 i/Q p [TS] . 

In the weak probe limit, the steady state solution of 
p 3 i of the A-type EIT system is given by [T5] 

P31 _ A p - A c + vy 

ftp ft2/ 2 -2(A p + ir/2)(Ap- A c + i 7 )' 1 ' 
which can be expanded as a power series of uj: 

g=£ &(i) ("-^ (5) 

where u p is the center frequency of the probe pulse. It 
has been shown that for an input Gaussian probe pulse 
with a width To, 

iW<) = ^oe-* 2/r », (6) 
the output probe pulse will take the form [TH] 

fW(t) = ft P o (j^j e-P-^ 2 /^, (7) 

where j3 = Lrjlm , = Lr\Re [b^'] and t w = 
\J t'$ + ALr/Im [&( 2 )] , corresponding to the attenuation 

constant, delay time and the broadened width of the out- 
put probe pulse, respectively. In deriving eq.Q, we re- 
tain the power series in eq. ([5| merely up to the quadratic 
terms and set A p = A c . Under the approximations of 
strong coupling field, ft 2 . S> 7A p and ft 2 ^ jT, M ' and 
&( 2 ) are purely imaginary numbers and b^ 1 ' is a real num- 
ber. The above conditions for achieving eq.Q are well 
justified in the EIT related experiments since the EIT 
bandwidth is proportional to ft;?/r and the frequency 
bandwidth of the probe pulse limits the maximum value 
of (uj — uj p ). Note that eq.([7]) is obtained based on the as- 
sumption that all atoms are kept fixed. Realistically, the 
atoms move freely with a velocity distribution obeying 
the Maxwcll-Boltzmann statistics at finite temperatures. 
To reconcile these two contradictory situations, eqs. ([!])- 
([2]) should be understood as the motion equations seen by 
any observer fixed on the moving atoms. To incorporate 
the Doppler shift into eq.Q in an uncomplicated way, we 
shall assume that the probe pulse propagates along the 
major-axis of the atomic cloud and the coupling field is 
applied with an angle 9 with respect to the major-axis. 
We also assume that the center frequencies of coupling 
and probe fields are on two-photon resonance in the lab- 
oratory frame, namely, A p = A c . Hence, for an atom 
moving with a velocity v = vi +v 2 in the non-relativistic 
limit, where v z is the velocity along the major-axis and 
vj_ is the velocity in the transverse direction, it expe- 
riences the Doppler shifts so that A p — > A p — k p v z and 
A c — > A c — (k c v z cos 9 + k c v± sin 9) [50] . As a result, the 
condition of two-photon resonance does not hold gener- 
ally, and the two-photon detuning, 8 = A p — A c , is max- 
imized when 9 = it. 



3 



In the presence of Doppler shifts, all fcW in eq. 
now turns out to be velocity-dependent, i.e., few — >. 
&^(v z , Statistically, can be replaced by an effec- 
tive coefficient obtained by taking the ensemble average 
of b( l '(v z ,v±) over the Maxwell-Boltzmann velocity dis- 
tribution at a given temperature T, 



g(0( v ) -> &« (T) 



1 



dv z dv ± U l) (v)e- v2/v % (8) 



where w s = y/2kBT/m is the one-dimensional root-mean- 
square velocity. It should be noted that the ensemble- 
averaging process in eq.Q is liable to the relaxation of 
P21 owing to the atomic thermal motion. To be consis- 
tent, in each (v), we let 7t = and keep 70 in 7 while 
evaluating the coefficients &W (T). Our effective the- 
ory is formulated based on the above ensemble-averaging 
process, in which the analytical forms of f2 c (T), fj(T) 
and 7 (T) are derived. The basic idea relies on the pre- 
sumption that, either the stationary medium character- 
ized by f2 c (T), fj (T) and 7 (T) or the Doppler-broadened 
medium characterized by the temperature-independent 
ft c , i] and 70, should lead to the same attenuation con- 
stant, delay time and the broadened width of the output 
probe pulse, namely, 



0_ = 
L 

r d (T) 
L 



6 (o) 



7,O c 



7o,fi c 



7,O c 



and 



rl{T) 



4L 



Ib^ 



7,n e 



(9) 



(10) 



(11) 



We evaluate eqs.(9| — ( 11 ) by making use of eqs.Q and 
together with the assumption A p = A c = for a 
stationary medium, and this yields 
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where h(v,0) — kv± sin# — kv z (1 — cos#), G(v z ) = 
-fcu, + iT/2 and 5 (v, 9) = 12 7o /i (v, 9) [ih (v, 0) - 7o ] + 
2z 7 of2;?. To further simplify the last three equations, 
we assume |sin0| -C 1 together with the assumptions 
ill > 7oI\ ill > 7 r > > 7o A P , r » 7 and 
k p m k c = k which are the typical conditions for the 
EIT experiments. The condition |sin0| <C 1 is generally 
true, since so far all EIT experiments have used either 
the nearly co-propagating (9 ~ 0) or counter-propagating 
(9 ss 7r) geometry. In our experiments, the central wave- 
length of the laser beam is 780nm, 70 ~ 5 x 10 _4 r, 
7 ~ 2 x 10~ 3 r, A p = 0, and 0.6r < fl c < 0.9T, 
where the spontaneous decay rate of the excited state 
is T = 2ir x 5.9 MHz. Obviously, the forementioned as- 
sumptions are valid and thus the effective parameters can 
now be solved as 
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and 



7(T) = 7o + 4/(0)§x 2 + 
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where / (9) = 1 - cos6> and x = sin (9/2) kv s T/V2il 2 c , 
which is the ratio of the effective Doppler width to that 
of the EIT window. In eq.( 17), the leading term of 7 (T) 
is the ground-state relaxation rate originating from the 
intrinsic effects and the second term is resulted from the 
Doppler shifts. When the probe and coupling fields are 
nearly co-propagating, i.e. 9 ~ 0, the Doppler effect is 
considerably suppressed, and this is why the EIT exper- 
iments at room temperatures are always studied in the 
co-propagating geometry |13j . On the contrary, when 
the probe and coupling fields are in the nearly counter- 
propagating geometry, i.e. 9 ~ n, the decoherence mech- 
anism is dominated by Doppler effects, and this is why 
the EIT-related experiments in the counter-propagating 
geometry are always studied in cold atoms [22] . 



III. RESULTS AND DISCUSSIONS 

To demonstrate the availability of the current scheme, 
_ v 2 , v we use the experimental data of SL as an input for 



[Ql/2-2[h(v,8) +i l0 ]G(v z )} 
fjirv 2 mg/2-47 3 + 27^ 
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cqs.( 15 )-( 17 ). To highlight the effect of thermal motion 



of atoms, we consider the case when the probe and cou- 
^4pling beams are arranged with 9 ~ 7r along the major 
axis of a cigar-shaped cloud of laser-cooled 87 Rb atoms 
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FIG. 2: (Color online) The squares and circles are the ex- 
perimentally measured intensities of the input and output 
probe pulses, respectively. The input probe pulses are scaled 
down by a factor of 0.1 in (a) and 0.01 in (b). The dashed- 
dotted lines represent the numerical simulations of the input 
probe pulses. The black lines represents the numerical re- 
sults of the output probe pulse by solving SL in the station- 
ary EIT medium using the effective parameters (f2 c , Lfj, 7) = 
(0.7T, 16.13r, 0.0143r) in (a) and (0.67r, 19. 981", 0.0227r) in 
(b). The gray dotted lines represent the best fitting by us- 
ing the scheme in [T7] with the parameters (T,Q C , Lrj) = 
(205/jA", 0.7ir, 16.5r) in (a) and {305fxK, 0.68I\ 20. 51") in (b). 
Clearly, the black and gray dotted lines completely merge in 
both (a) and (b). 



[2T] . The experimental values of the coupling Rabi fre- 
quency, opacity and relaxation rate of the ground-state 
coherence (denotes as CI* , Lrf and 7q ) can be estimated 
by the method described in [22] . We apply the numer- 
ical scheme of p~7| , which considers all velocity groups 
of the atoms in the calculations, to determine the tem- 
perature T of the atomic medium and get the best fit- 
ting of fi c and Lrj as the gray dotted lines plotted in 
Fig. 2(a) and (b). To be generic, we consider two inde- 
pendent sets of experimental data of SL with the esti- 
mated values (Q*,L?7*,75) = (0.687r, 15r, 0.0005r) and 
(0.69r, 20.5r, 0.0005r). The best fitting comes up corre- 
spondingly with (T, O c , Lrj) = (205/xK, 0.7LT, 16.5r) and 
(305/xfT,0.68r,20.5r). Inserting the above best fitted 



parameters into eqs.(15)-(17l, the effective parameters 
are found to be (fi c , Lfj, 7) = (0.7r, 16.13r, 0.0143r) and 
(0.67r, 19.98r,0.0227r). Using fl c) Lfj and 7 in eqs.Q- 
([3]) , we numerically calculate the output probe pulse as a 
function of time and find it is in good agreement with the 
experimental data as the black lines plotted in Fig. 2(a) 
and (b). The validity of eqs.(15)-(17) is hence demon- 
strated. 
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FIG. 3: (Color online) The fractional changes, SQ C = 
(O c — Q c ) /f2 c and 5fj = (77 — fj) /rj are plotted in (a) and 
(b) as functions of temperatures and £l c with 9 = tt. 



The proximity of Q c and fj to their zero-temperature 
counterparts can be examined by defining the fractional 
changes, SCl c = (fl c — J7 C ) /fi c and Sfj — (77 — fj) /rj, 
which are shown in Fig. 3 as functions of fl c and T within 
the experimentally accessible regime. Obviously, the 
small values of 5fl c and Sfj shown in Fig. 3 indicates that, 
Q c « Q c and r\ fj, a conclusion that is consistent with 



the perturbative results in eqs.( 15 1-( 16 ). We thus expect 
that the atomic thermal motion affects the ground-state 
relaxation rate mostly. In Fig. 4, we compare the mea- 
sured ground-state relaxation rate with the theoretical 
prediction, eq.(17). Here the measured relaxation rate 



is determined by fitting the experimental data with the 
numerical results from eqs.(l)-(3). In Fig. 4, the theoret- 
ical predictions agree with the measured data very well 
and we see that the averaged temperatures of our exper- 
imental setup of Lrf = 15r and 20. 5r are found to be 
T = 200/jK and 330/jK, respectively. 



After demonstrating the validity of eq.(17l, we study 
how the decoherence rate 7 varies with f2 c , T and 9. 
Giving 9 = tt, we first depict 7 as a function of f2 c and 
T in Fig. 5. Here we extend the ranges of Cl c and T so 
that the conditions of stationary light pulses and those 
using the counter-propagating beam geometry with laser- 
cooled Rb atoms can be included. We see that 7 decays 
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FIG. 4: (Color online) The comparisons between the mea- 
sured 7 and the theoretically predicted 7 for Lr;* = 20. 5r, 
15F in (a) and (b), respectively. The circles are the mea- 
sured data and the solid lines are the theoretical predictions 
at T = 330^iK, 200/xK in (a) and (b). The corresponding 
temperatures of the measured data in (a) range from 290/^K 
to 370/^K while those in (b) range from 175/iK to 230/iK. 



with increasing intensity of the coupling field as shown in 
Fig. 5. This is because the width of the transparency win- 
dow depends on the intensity of the coupling field, and an 
increasing fl c can lower the absorption of the probe pulse, 
so that 7 would become smaller. The influence of the ori- 
entation of the coupling beam and the effect of Doppler 
width on 7 is investigated by plotting 7 as a function of 
9 and k p v s /Q, 2 in Figs. 6-7, which were obtained based on 
the numerical calculations of eqs.(12l — (14). The condi- 



tions of most experiments utilizing small separation angle 
between two propagating beams in lase-cooled atoms are 
met in the ranges, 1° < 6 < 21°, O.Oir -1 < k p v s /Q 2 c < 
r _1 , as shown in Fig. 6. On the other hand, those ex- 
periments utilizing the nearly co-propagating beam ge- 
ometry for hot atoms are mostly achieved in the ranges 
0.1° < 6 < 2.1°, 10T- 1 < kpV s /Q 2 c < lOOOr- 1 , as shown 
in Fig. 7. Taking our experimental conditions ( 87 Rb 
atoms and A = 780nm) as an example, if T = 100/iK 
then we have k p v s /i} 2 ss 0.03F" 1 . Figures 5-7 thus pro- 
vide useful information to determine the relaxation rate 
for those EIT-related experiments in which the Doppler- 
broadening cannot be ignored. 



IV. CONCLUDING REMARKS 

In conclusion, we have theoretically studied the effects 
of atomic thermal motion on the dynamics of SL in a 
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FIG. 5: (Color online) Contour plot of 7 with 9 — ir. The 
vertical axis ranges from T — 1/jK to 1000/iK and the hori- 
zontal axis ranges from fi c = 0.4r to 1.2r. The values of the 
contour lines are from 0.12 (the upper left), 0.1, 0.08, 0.06, 
0.04, 0.02, 0.01,0.005 (the lower right). 
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FIG. 6: (Color online) The plot of 7 as a function of and 
kpV s /Q^, in logarithm scale. The vertical axis ranges from 
8 — 1° to 21° and the horizontal axis ranges from k p v B /Q^. 
( ) J ) 1 f to r _1 . The values of the contour lines is from 10 -3 
(the upper right) to 10" 10 (the lower left) with a spacing of 
one order of magnitude. 



Dopplcr-broadencd EIT medium. By extending the re- 
sults in [TH] , we have obtained the attenuation constant, 
the delay time and the broadened width of the output 
probe field at finite temperatures. We also have derived 
a set of parameters, f2 c , fj and 7, which are temperature- 
dependent that serve as effective thermal sources in the 
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optical Bloch equation for a stationary EIT medium. Our 
approach provides an approximate and convenient way to 
determine the relaxation rate of SL instead of appealing 
to the more complicated multiple-velocity calculations. 
This effective theory is not only valid for the cold atomic 
medium but also applicable to the atomic medium at 
room temperatures. 
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FIG. 7: (Color online) The plot of 7 as a function of 9 and 
kpV B /Qc. The vertical axis ranges from 6 — 0.1° to 2.1° and 
the horizontal axis ranges from from k p v s /0,^. = 10r _1 to 
1000r _1 . The values of the contour lines are from 10° (the 
upper right) to 10 -8 (the lower left) with a spacing of one 
order of magnitude. 
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